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COUPLED  AMPLIFIER  MODULE  FEED  NETWORKS 
FOR  PHASED  ARRAY  ANTENNAS 


INTRODUCTION  AND  SUMMARY  OF  RESULTS 

The  high  cost  of  phased-array  antennas  prevents  their  use  in  many  applications  where  electronic 
scanning  could  provide  large  performance  improvements  relative  to  mechanically  scanned  antennas. 
One  reason  for  the  high  cost  of  conventional  phased-arrays  is  their  parallel  construction:  each  radiating 
element  has  its  own  phase  control.  Previous  approaches  to  reducing  system  cost  by  even  modest  reduc¬ 
tions  in  the  number  of  phase-shifters  result  in  objectionable  increases  in  sidelobe  levels  [1].  Although 
the  feed  networks  discussed  in  this  report  allow  for  greatly  reduced  numbers  of  phase-shifters,  there  are 
other  components  that  must  be  added.  Thus  the  cost-effectiveness  of  the  new  approach  is  an  issue  for 
further  study. 

Figure  1  depicts  a  linear  array  corporate  feed,  illustrating  the  one  phase-shifter/element  require¬ 
ment  that  characterizes  conventional  design  [2].  Figure  2  shows  the  functionally  equivalent  coupled 
amplifier  module  (CAM)  feed.  Phase  control  is  provided  only  at  the  edges  of  the  CAM  feed,  and 
requires  just  one  phase  shifter.  As  shown  in  Fig.  3,  each  module  in  the  CAM  feed  consists  of  a  hard- 
limited  amplifier  whose  input  is  the  sum  of  sampled  outputs  from  the  two  adjacent  modules. 


RADIATING 

ELEMENTS 


PHASE  SHFTERS 


HYBRID  JUNCTION 
CORPORATE  FEED 


Fig.  2  —  Four-element  linear  array  employing  coupled  amplifier  module  feed 
Manuscript  approved  November  8,  1982. 
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Fig.  3  —  Amplifier/limiter  module  for 
coupled  amplifier  module  feed 


Since  (he  outputs  of  the  amplifier/limiter  modules  are  all  of  unit  amplitude,  all  information  is 
contained  in  their  phases,  0„(t),  where  n  is  the  module  index  and  r  is  the  time.  It  is  shown  that  the 
dynamics  of  the  CAM  network  are  governed  by  a  simple  diffusion  equation  on  the  phases  0„(t). 

The  main  body  of  this  report  is  devoted  to  developing  solutions  to  the  phase  diffusion  equation. 
Consideration  is  given  to: 

•  steady-state  phase  distributions  for  fabricationally  perfect  modules,  i.e.,  modules  for  which  the 
power  dividers  effect  exactly  equal  two-way  power  division; 

•  the  complete  transient  phase  distribution  for  perfect  modules;  and 

•  steady-state  phase  distributions  for  fabricationally  imperfect  modules,  including  those  for  which 
the  imperfection  can  only  be  characterized  statistically. 

The  phase  diffusion  equation  can  be  solved  exactly  when  the  modules  are  fabricationally  perfect, 
with  the  following  results. 

•  While  the  CAM  network  establishes  a  linear  phase  progression  in  the  steady  state,  the  phase 
difference  between  the  outputs  of  adjacent  modules  is  limited  to  the  interval  [0°,  90°).  If  the 
phase  gradient  were  ever  to  reach  90°  per  module,  the  signals  summed  at  the  input  to  the  hard- 
limited  amplifier  would  be  1 80°  out  of  phase  and  therefore  would  add  to  zero.  Since  it  is  thus 
impossible  to  develop  phase  gradients  greater  than  ±90°  per  module,  two  modules/element  are 
required  to  achieve  full  hemispheric  scanning  capability  (cf.  Fig.  2).  (Only  one  module  per  ele¬ 
ment  may  be  required  for  limited  scan  applications,  however.) 

•  The  exact  transient  response  takes  (he  form  of  a  modal  expansion  in  which  each  of  the  spatial 
modes  has  a  characteristic  decay  time.  Since  the  transient  response  is  dominated  by  the  spatial 
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mode  having  the  slowest  decay  time,  a  simple  and  highly  accurate  one-mode-approximation 
results. 

•  The  settling  time  of  the  CAM  network  increases  as  the  square  of  the  number  of  modules.  Thus, 
satisfaction  of  bandwidth  requirements  may  necessitate  introducing  additional  "speed-up"  phase 
shifters  spaced  at  regular  intervals  along  the  feed  network. 

The  analysis  for  fabricationally  imperfect  modules  permits  only  approximate  solutions  for  the 
steady-state  phase  distributions.  Defining  the  phase  gradient  error  as  the  difference  between  the  actual 
phase  gradient  and  the  ideal  constant  gradient,  an  expression  is  developed  for  the  phase  gradient  error 
for  arbitrarily  nonuniform  CAM  structures.  Assuming  that  the  nonuniformities  from  one  module  to 
the  next  are  independent,  equal-variance  random  variables,  an  expression  for  the  root-mean-square 
(RMS)  phase  gradient  error  is  then  derived.  It  follows  from  this  analysis  that  the  RMS  phase  gradient 
error  increases  as  the  square  root  of  the  number  of  modules  in  the  CAM  network;  also,  the  error  is 
proportional  to  tan  A0,  where  A<£  is  the  nominal  phase  increment  per  module.  An  example  is  given 
for  a  9-module  feed  for  which  a  1%  RMS  error  in  the  power  splits  gives  rise  to  about  a  1°  RMS  phase 
gradient  error,  for  an  assumed  nominal  phase  gradient  of  45°  per  module. 

Finally,  it  is  shown  that  the  CAM  technique  can  be  applied  to  2-D  antenna  arrays.  The  CAM 
feed  network  for  a  2-D  rectangular  array  of  arbitrary  size  requires  just  three  phase-shifters. 

In  most  practical  antennas  the  amplifier/limiter  modules  feeding  the  radiating  elements  would  be 
followed  by  power  RF  amplifiers.  In  receive  mode,  the  coupled  module  structure  would  be  used  as  a 
local  oscillator  chain  to  convert  the  antenna  outputs  to  intermediate  frequency  (IF). 

MODEL  FORMULATION 

Figure  4  shows  an  idealized  model  for  the  individual  amplifier/limiter  modules.  The  principal 
idealization  is  in  representing  the  hard  limiter  as  a  time-variable  amplifier  whose  gain  /V„(t)  automati¬ 
cally  adjusts  itself  to  provide  a  unit  amplitude  output,  i.e., 

x„(t)  -  exp  j0„(r).  (1) 

where  n  is  the  module  index,  and  the  normalized  time  r  is  defined  as 

t  -  (r/D, 

where  T  is  the  delay  time  per  module  (cf.  Fig.  4).  The  gain  N„( t)  is  a  function  of  time  and  is  generally 
different  from  one  module  to  the  next, 

N„(t)  &  Nm(r),  rt  ^  m.  (2) 

We  preserve  sinusoidal  waveforms  by  implicitly  assuming  that  the  hard  limiter  in  Fig.  3  is  followed  by  a 
high-Q  tuned  circuit. 

The  gains  ,  and  F„+\  in  Fig.  4  are  determined  from  the  power  splits  effected  by  the  two  power 
dividers  in  Fig  3  and  are  respectively  referred  to  as  the  backward-  and  forward-coupling  coefficients. 
We  generally  design  the  power  splitters  to  achieve  the  ideal, 

-  F„.  (3) 

Departures  from  the  ideal  are  unfortunately  inevitable  due  to  fabricational  imperfections  in  the  power 
splitters. 

From  Figs.  2  and  4, 

x„(r  +  1)  -  N„(r)lB„xn+,(r)  +  F„x„.,( r)). 
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Fig.  4  —  Idealized  model  of  amplifier/limiier 
module  depicted  in  Fig.  3 


Thus,  from  Eq.  (1), 

*"'(T+"  -  Njr)lBHei0^{T)  +  Fnej9-'M) 

where  n  =  1,  2,  ...  (M  —  1).  Solution  of  Eq.  (4)  generally  requires  the  specification  of  boundary 
phases  at  n  —  0  and  n  =  A/, 


(4) 


and  the  initial  phases 


90(t)  -  0 

(5a) 

1 »u 

-  —  r 

•  2  n. 

T  <  0 

-  |,t 

~  A/A<£  —  s 

■  2ir, 

T  ^  0, 

(5b) 

- 

»-  1,  2,  . 

..  A/, 

r  <  0. 

(6) 

The  terms  r  ■  2n  and  s  •  2ir  appearing  in  Eq.  (5b),  where  r  and  s  are  integers,  reflect  the  fact  that  the 
phase  difference  (0M  -  0O)  established  by  the  phase-shifter  (cf.  Fig.  2)  is  physically  determinate  only  to 
within  an  integral  multiple  of  27r.  Also,  we  note  that  the  number  of  modules  in  the  feed  network  is 
equal  to  (M  -  1).  (For  example,  M  —  6  for  the  network  depicted  in  Fig.  2.) 


For  future  use  we  define  the  quantity  aM(0)  as  the  step  change  in  boundary  phasing  imposed  by 
the  phase-shifter  at  t  -  0, 

ctjVf(O)  “  modjjffflj^’  —  &i{)  ”  A/(A<£  —  A<£q).  (7) 

The  function  modlr  (•)  appearing  in  Eq.  (7)  expresses  its  argument  in  modulo  -  2 ir,  i.e.,  as  a  number 
between  0  and  2ir. 
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THE  PHASE  DIFFUSION  EQUATION 

It  is  shown  in  Appendix  A  that  Eq.  (4)  can  be  recast  as 


9„( t  +  I) 


1 


0„+i (r)  +  6>n_,(r)  -  e„  tan  y  [0„+,( r)  -  0„_ ,(t)]  , 


(8) 


which  we  refer  to  as  the  phase  diffusion  equation.  The  parameter  e„  provides  a  measure  of  fabrica¬ 
tion!  imperfection,  being  defined  as 


Fn-B„ 

Fn  +  BA 


(9) 


Equation  (8)  is  derived  in  Appendix  A  subject  to  the  assumption  that 

y  € „  tan  y  10„+|(t)  -  «  1.  (10) 

Thus,  Eq.  (8)  is  only  approximately  correct  when  ^  0;  however,  the  simpler,  less  general  equation 

that  results  when  e„  —  0  is  exactly  equivalent  to  Eq.  (4)  when  F„  —  B„. 

EQUILIBRIUM  PHASE  DISTRIBUTIONS 

Analysis 

Equilibrium  solutions  of  Eq.  (8)  are  obtained  by  recognizing  that 

lim  9„(t  +  1)  «  lim  9„(r)  S  9n.  (11) 

the  final  equality  being  a  convenient  notational  convention.  With  the  additional  definition 

(12) 

it  follows  from  Eqs.  (8)  and  (11)  that  in  the  infinite-time  limit 

A 0„+,  -  AO,  +  €„  tan  y  (A0„+l  4-  A 9„).  (13) 

If  we  are  to  have  a  uniform  phase  gradient, 

A0,  -  A0„+|  —  constant,  (14) 

it  follows  from  Eqs.  (9)  and  (13)  that  we  must  have  equal  forward  and  backward  coupling, 

Fn  “  Bn,  (15) 

i.e.,  «„  —  0.  Feed  networks  for  which  Eq.  (15)  obtains  are  referred  to  as  "fabricationally  perfect."  If  F„ 
&  B„,  i.e.,  if  «„  0,  for  one  or  more  values  of  n,  the  network  is  referred  to  as  "fabricationally  imper¬ 

fect." 


From  Eqs.  (5)  and  (12), 


and 


i  A0m  -  (•„  -  90)  -  9„, 

/»•  1 

£  A0m  -  (0W  -  90)  -  M£nt>  -  s  ■  2n, 


(16) 


(17) 
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where  s  is  an  integer.  The  term  s  •  2tt  appearing  in  Eq.  (17)  reflects  the  Tact  that  the  phase  difference 
(9 \i  -  90 )  established  by  the  phase  shifter  (cf.  Fig.  2)  is  physically  determinate  only  to  within  an 
integer  multiple  of  2-ir. 

From  Eq.  (14), 

£  \8m-n\e„,  n  -  I,  2.  ...  M.  (18) 

m—  1 


The  equilibrium  phase  distribution  for  fabricationaily  perfect  CAM  networks  follows  from  Eqs. 
(16),  (17),  and  (18). 


A0„(°o)  —  A</>  -  (s  ■  lir/M) 

n\e„(°o)  n  -  2,  ...  M. 


(19) 


where  s  is  an  integer. 


To  complete  the  specification  of  the  physical  problem,  and  to  provide  necessary  information  for 
calculating  the  indeterminate  integer  s  appearing  in  Eq.  (19),  we  must  specify  values  for  the  initial 
phase  distribution,  Eq.  (6).  Assuming  that  the  network  is  initially  in  an  equilibrium  state  we  may  write 

6»b(t)  -  n  -  1.  2,  ...  M,  r  <  0.  (20) 


For  initial  conditions  and  boundary  conditions  specified  by  Eqs.  (20)  and  (5),  respectively,  it  has 
been  found  by  computer  simulation  that  the  only  values  of  sever  to  occur  are  s  “  0  and  s  —  1. 

Figure  5  shows  the  initial  phase  distribution,  Eq.  (20),  and  the  two  possible  final  phase  distribu¬ 
tions  (Eq.  (19)  with  s  =-  0  and  s  =  1).  We  now  summarize  some  results  from  Appendix  B  addressing 
the  question  of  which  of  these  two  final  phase  distributions  is  ultimately  established. 


For  positively  coupled  lattices,  F  >  0, 


Jo,  0  <  aM(0)  <  (w  -  2A$0) 
s  “  ( 1 ,  otherwise. 

For  negatively  coupled  modules,  F  <  0, 

0,  0  ^  GfjifC 0)  K.  (3ir  2A^0) 
s  **  |l ,  otherwise. 


(21) 


(22) 


We  recall  that  the  quantity  aM(0)  appearing  in  Eqs.  (21)  and  (22)  was  defined  in  Eq.  (7)  as  the  step 
change  in  boundary  phasing  imposed  by  the  phase-shifter  at  time  r  —  0. 

To  summarize,  the  equilibrium  phase  distribution  9„  (or  A0„)  is  calculated  subject  to  two  assump¬ 
tions: 

•  We  assume  that  the  network  is  initially  in  an  equilibrium  state  characterized  by  the  initial  phase 
gradient  A 9„(t)  —  A<£o.  r  <  0. 

•  We  assume  that  the  setting  of  the  phase  control  in  Fig.  2  is  changed  by  an  amount 
mod2„(04  “  9'm)  =  ««(0),  having  a  value  between  0  and  2 ir. 
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Fig.  5  —  Phase  distributions  for  an  8-module  linear  lattice  ( M  —  9). 
The  initial  phase  distribution  0'  has  a  phase  gradient  of  A^0  per 
module.  The  initial  distribution  evolves  into  one  of  two  possible  final 
distributions,  «/,  referred  to  as  the  s  —  0  and  s  *•  1  equilibrium 
phase  distributions.  The  s  “  0  distribution  is  established  if  0  < 
modj,  (»{,  -  fiji,)  <  (»  -  2A40) 


The  network  will  settle  into  a  new  equilibrium  phase  distribution  that  can  be  calculated  as  follows. 

•  The  value  of  sis  computed  from  Eq.  (21)  if  F >  0  or  from  Eq.  (22)  if  F  <  0. 

•  The  quantity  A<£  is  computed  from  the  equation 

A«*  -  A0O  +  l a„(0)/M).  (23) 

Since  o*/(0)  is  always  positive  by  definition  (cf.  Eq.  (7)),  Eq.  (23)  ensures  that  A <t>  >  A^0. 
always. 

•  The  values  for  b<f>  and  sare  substituted  into  Eq.  (19)  to  obtain  the  new  equilibrium  phase  distri¬ 
bution. 

Appendix  B  provides  expanded  discussion  of  the  equilibrium  phase  distributions. 


Phasing  Limitations  and  an  Alternative  Module  Structure 

Inspection  of  Fig.  $  shows  that  development  of  the  s  —  0  phase  distribution  corresponds  to  the 
evolution  of  a  larger  phase  gradient  than  existed  initially, 

A0„  >  A*0.  s  -  0.  (24) 
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It  follows  from  Eqs.  (21)  and  (24)  that 

aM(0)  <  (n  —  2A<£o)>  (25) 

from  which,  together  with  Eqs.  (19)  and  (7),  we  can  show  that 

(A0„  —  A<^0)  <  (w  —  2A£o)/ M-  (26) 

Equation  (26)  indicates  that  ihe  achievable  increment  in  phase  gradient  becomes  ever  smaller  as  the 
initial  gradient  approaches  90°  per  module.  Thus,  phase  gradients  |A6»J  >  90°  per  module  are  unattain¬ 
able  with  positively  coupled  lattices  (F  >  0).  Similarly,  we  find  that  phase  gradients  0  <  |A0„|  <  90°  per 
module  are  unattainable  for  CAM  lattices  with  negative  coupling  ( F  <  0).  Appendix  B  provides  a 
more  detailed  discussion  of  the  phasing  limitations. 

One  method  of  achieving  full  hemispheric  scanning  capability  (up  to  ±180°  phase  differential 
between  antenna  elements)  is  to  allocate  two  modules  per  radiating  element  (cf.  Fig.  2).  Alternatively, 
only  one  module  per  antenna  need  be  used  if  the  module  is  designed  as  shown  in  Fig.  6.  As  compared 
with  the  original  Fig.  3  module,  the  Fig.  6  module  incorporates  a  180°  1-bit  phase-shifter.  The 
modified  CAM  network  is  made  to  simulate  the  original  network  for  small  phase  gradients,  0  ^  A0  < 
90°  per  module,  by  setting  all  I -bit  phase-shifters  to  their  0°-state.  Scanning  at  large  angles,  requiring 
90°  <  A 0„  <  180°  per  module,  is  achieved  by  setting  all  1-bit  phase-shifters  to  their  180°-state 
(corresponding  to  negative  coupling,  F  <  0).  Thus,  utilizing  the  Fig.  6  module  rather  than  the  Fig.  3 
module  we  may  replace  all  nonradiating  terminations  in  Fig.  2  with  radiating  elements  and  still  retain 
full  hemispheric  scan.  Of  course,  the  expense  saved  in  halving  the  number  of  modules  may  be  more 
than  completely  offset  by  the  increased  cost  per  module. 


Fig.  6  —  Amplifier/limiter  module  incorporating  I -bit  180°  phase-shifter, 
permitting  full  hemispheric  scan  with  just  one  module  per  radiating  ele¬ 
ment 
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TRANSIENT  PHASE  DISTRIBUTION  AND  SETTLING  TIME 


We  presently  assume  that  the  CAM  lattice  is  fabricationally  perfect,  i.e.,  that  Eq.  (15)  obtains  for 
all  values  of  n.  It  follows  from  Eqs.  (8),  (9),  and  (IS)  that 

0„(t  +  1)  -  y  (0„+1(t)  +  0„_|(t)].  (27) 

Equation  (27)  is  solved  in  Appendix  C  by  the  separation  of  variables  method.  The  exact  solution  is 

Mr)  -  M°°)  +  2UM“>  -  A*0]  1 if  (-°*[ 'shi {(kn/Mi  fc°S  (kn/M)]T+ (W 

where  [M/2]  denotes  the  integer  part  of  (A//2),  quantity  A0„(°°)  is  the  equilibrium  phase  distribution, 
Eq.  (19),  and 


M  =  y  ll  +  (-1)"+T+W) 


0  ,  (n  +  r  4-  M)  odd 
ll  ,  (/»  +  r  4-  M)  even. 


(29) 


The  factor  [A0„(<»)  -  A<£0)  appearing  in  Eq.  (28)  may  be  calculated  from  Eqs.  (19)  and  (23)  as 

IA<U°°)  ~  A*o)  -  (aw(0)  -  s2rr]/M  (30) 

where  aw(0)  is  the  step-change  in  edge  phasing,  as  defined  by  Eq.  (7),  and  s  is  given  by  Eq.  (21)  if 
F>  Oor  Eq.  (22)  if  F<  0. 

Values  of  0„(r)  calculated  from  Eq.  (28)  have  been  compared  with  the  results  of  iterative  numer¬ 
ical  solutions  of  Eq.  (4).  The  exact  equality  of  0„  values  calculated  by  these  two  totally  different 
methods  is  taken  as  a  validation  of  Eq.  (28). 

In  the  limit  ( kn/ M )  <<  1,  valid  for  small  fc, 

Icos  ( kir/M)Y  =s  e  t/t‘,  (31) 

where  the  time  constants  rk  are  given  by 

r*  =  |  ( M/k )2.  (32) 

It  follows  from  Eqs.  (31)  and  (32)  that  the  k  ^  2  terms  in  Eq.  (28)  decay  much  more  rapidly  with  time 
than  does  the  k  =•  1  term.  Thus,  we  expect  the  long-time  behavior  of  the  phase  distribution  to  be 
dominated  by  the  A  =  1  term, 

Mr)  -  0„(°o)  -  2[AMoo)  -  A<ft0I  ~'n  {f!nl^\e~,r^)h5'  (33) 

[sin  (ir/M) 

where  n  is  again  given  by  Eq.  (29),  and  the  settling  time  rs  is  defined  as 

Tj  =  T]  —  M2/S.  (34) 

A  simple  single-mode  approximation  also  may  be  obtained  for  the  phase  gradient  A 9n,  Eq.  (12), 

A 9„(t)  —  A0„(°°)  —  2[A0„(<»)  —  A<£ol  cos  [(n  +m  -  \)n/M]e  t/t'.  (35) 

The  single  mode  approximations,  Eq.  (33)  and  (35),  are  highly  accurate  for 

T  >  rs. 

We  have  compared  the  approximate  solution,  Eq.  (35),  with  the  complete  analytic  solution 
obtained  from  Eq.  (28)  and  with  iterative  computer  solutions  of  Eq.  (4).  The  sample  calculations 
presented  below  assume  that 
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(A/,  \<t> o,  a  ,/ (0) )  =  (9,  60°,  45°).  (36) 

Thai  is.  we  assume  an  8-module  feed  initially  in  an  equilibrium  slate  described  by  Eq.  (20)  with  a  uni¬ 
form  phase  gradient  of  60°  per  module.  We  further  assume  that  the  equilibrium  prevailing  for  r  <  0  is 
disturbed  at  r  =  0  by  step-increasing  t-he  phase  control  in  Fig.  2  by  a^lO)  =  45°.  We  find  from  Eqs. 
(21)  and  (23)  that  the  network  will  eventually  settle  into  a  new  equilibrium  phase  distribution  given  by 
Eq.  (19)  with 

A0„(«)  »  A<ji  =  65°.  (37) 

For  this  calculation  (and  all  other  calculations)  the  numerical  solution  was  in  exact  agreement 
with  Eq.  (28). 

Figure  7  shows  the  phase  distribution  as  a  function  of  module  index  n  at  two  values  of  time 
(r  »  32  and  t  =  48),  for  the  specified  parameter  values,  Eq.  (36).  Approximation  A0„  is  indistinguish¬ 
able  from  the  exact  solution  on  the  scale  of  Fig.  7. 


Fig.  7  —  Phase  distribution  for  an  8-module  feed  at  two 
values  of  time  (r  -  22  and  t—  48).  Assumed  parameter 
values  were  (M,  A(p0.  «w(0))  •>  (9.  60°,  45 °).  In  the 
infinite-time  limit  a  uniform  phase  gradient  is  established, 
40„(«o)  -  -  65- 


Figure  8  depicts  the  transient  development  of  A0i(t)  obtained  by  the  exact  methods  of  solution. 
Figure  9  shows  the  difference  between  the  exact  and  approximate  solutions  for  A0|(r)  as  a  function  of 
time.  We  see  from  Fig.  9  that  the  error  incurred  in  using  Eq.  (35),  rather  than  the  complete  solution, 
is  less  than  0.13°  for  r  >  rs  *  16. 

MODULES  WITH  FABRIC  ATIONAL  IMPERFECTIONS 

Approximate  Formulation  for  the  Equilibrium  Phase  Distribution 

Up  to  this  point  in  our  analysis  we  assume  that  all  components  in  the  CAM  network.  Fig.  2,  can 
be  fabricated  to  arbitrarily  precise  specifications.  In  particular,  we  assume  that  the  power  splitters  have 
been  fabricated  to  achieve  perfectly  equal  forward  and  backward  coupling,  Eq.  (15).  It  is  of  some 
interest  to  quantify  how  component  imperfections  degrade  the  phasing  accuracy  of  the  CAM  network. 
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exact  methods  of  solution:  the  complete  analytical  solution,  r  >  r,. 

Eq.  (28),  or  by  numerical  solution  of  Eq.  (4). 


We  presently  relax  Eq.  (15)  somewhat,  assuming  instead  that 

0  <  |«J  «  l,  (38) 

for  at  least  one  value  of  n  (cf.  Eq.  (9)).  Our  interest  is  in  obtaining  the  solution  of  Eq.  (13).  We  have 
shown  previously  that  when  e„  ”  0,  n  “  1,  2,  ...  (A/  -  1),  the  equilibrium  phase  distribution  is  a 
uniform  phase  gradient,  as  desired, 

A0°(°o)  -  A<*o  +  [aw(0)  -  s  •  2ir]/M,  (39) 

or,  equivalently. 


A0«(oo)  -  A*  -  s(2ir/M).  (40) 

The  superscript  0  on  A0®(<»)  serves  as  a  reminder  that  Eqs.  (39)  and  (40)  are  derived  subject  to  the 
assumption  that  t„  —  0.  We  now  define  a  new  quantity  A<£„  that  measures  the  degree  to  which  the 
actual  phase  gradient  A 9„  departs  from  the  ideal,  i.e., 

Atfa  -  A0„(oo)  “  A0“(~),  (41) 

with  A0®(<»)  given  by  Eq.  (39). 
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Appendix  D  shows  [hat  Eq.  (13)  has  ihe  approximate  solution 

I  i.*p-  I4/-  n  =  1,  2.  ...  M.  (42) 

(  M  r- 1  p-r  p-n 

In  the  remainder  of  this  section  we  assume  that  s  =  0,  allowing  us  to  replace  A0®(«)  by  A <t>  in  Eq. 
(42). 

Feeds  with  Identical  Imperfect  Modules 

If  all  the  modules  in  the  CAM  network  are  identical,  even  though  they  are  imperfect,  we  have 

(43) 

independent  of  the  module  index  n.  Subject  to  Eq.  (43),  Eq.  (42)  becomes  (cf.  Eq.  (D26)) 

A <t>„  -  e  [n  -  ( M  +  l  )/2}  tan  A^.  n  —  1 ,  2,  ...  M,  (44) 

wh“£  «  -  2C F-  B)HF  +  SI. 

In  Fig.  10  we  present  values  predicted  by  a  computer  solution  of  Eq.  (4),  for  the  parameter  values 

(«,  A*.  M)  -  (Hr2,  65 9).  (45) 

We  observe  from  Fig.  1 1  that  Eq.  (44)  agrees  fairly  well  with  the  numerical  solution.  In  Appendix  D 
we  calculate  a  second-order  correction  to  Eq.  (44)  that  results  in  an  even  closer  fit  to  the  numerical 
solution. 


Fig.  10  —  Equilibrium  phase  error  for  an  8-module  feed  with  idemical 
imperfect  modules,  obtained  by  numerical  solution  of  Eq.  (4).  Parameter 
values  assumed  are  (*.  A <t>,  M)  “  (0.01,  65°,  9). 
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(I) 


Fig.  II  —  Difference  between  exact  value  of  A$„,  as  given  in  Fig.  10,  and 
value  of  A$„  obtained  by  two  approximate  formulations.  The  first-order 
approximation  A^„  is  given  by  Eq.  (44);  the  second-order  approximation 
is  given  by  Eq.  (D41). 


We  note  from  Fig.  10  that  the  phase  increment  per  module  A0„  varies  from  a  minimum  of 

A#  |  =  (Atf,  +  A  <{>)  =  60  ° 

to  a  maximum  of 

A0/n  “  (A <j>M  +  A<£)  —  70°, 

across  the  network.  Thus,  a  systematic  imperfection  of  1%  in  the  power  split  ratio  gives  rise  to  a  phase 
inaccuracy  of  ±5°.  However,  smaller  phase  errors  are  obtained  for  smaller  phase  gradients.  A <t>  <  65°, 
due  to  the  tan  A 4>  factor  in  Eq.  (44). 


Defining  a  mean-square  phase  error  as 

*  T7  Z  2 ■ 

M  n-\ 

it  follows  from  Eq.  (44)  that  (in  radians) 

f 

1 


rms  -  «  tan  Ad> 


12 


1/2 


M. 


For  our  previous  example,  Eq.  (45),  we  calculate 

A*rim-3.2°. 


(46) 


(47) 

(48) 


Random  Imperfections 

Rather  than  take  as  a  known  function  of  n  (e  g.,  Eq.  (43))  we  now  assume  that  the  coupling 
coefficients  F„  and  B„  are  independent,  identically  distributed  random  variables.  The  quantities  «„  are 
thus  also  independent  and  identically  distributed, 

w 

£{«,«*!  -  j0 

where  £{•}  is  the  statistical  expectation  operator. 
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P  -  k 
p  ^  k' 


(49) 
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We  show  in  Appendix  E  that  the  following  approximate  relationship  obtains  oeiween  the  statistics 
of  «  and  those  of  F. 

a,  =  V2  x  HT3  A.  (50) 

where  A  is  the  percentage  fabricational  error, 

A  =  \00{(Tt/nf).  (51) 


The  quantity  <rt  in  Eq.  (50)  is  the  standard  deviation  of  e„;  quantities  07  and  in  Eq.  (51)  are  the 
standard  deviation  of  F„  and  the  mean  value  of  F„,  respectively. 


The  variance  of  the  phase  error  is  defined  as 

<r2  =  F(A<£2)  -  (FTA^,,)]2. 


(52) 


It  is  shown  in  Appendix  D  from  Eqs.  (42),  (49),  and  (52)  that  (in  radians) 


tr„  =  <re  tan  &<t>  |-jy  1  +  l2AT2lrt  -  (M  +  l)/2]2}l/2. 


li/J 


(53) 


Equation  (53)  displays  the  expected  symmetry  about  n  =  ( M  +  l)/2.  The  minimum  value  of  cr„  is 
obtained  at  the  center  of  the  array;  from  Eq.  (53) 

1  41/2 

min  ar„  =  at  tan  A^l 


1 


12 


M1/2. 


(54) 


Equation  (54)  bears  an  interesting  resemblance  to  Eq.  (47);  however,  it  should  be  remembered  that 
■  was  defined  on  the  basis  of  a  spatial  average,  Eq.  (46),  while  the  values  of  cr„,  n  =  1,  2,  ...  M, 
are  established  as  ensemble  coverages  over  the  statistics  of  F„. 


From  Eqs.  (50),  (51)  and  (53)  we  obtain  (in  degrees) 


<t „  -  0.234  A  tan  A0  M,/2 


I  +  12AT2{n  -  (M  +  l)/2]2  ,  «  —  1,  2,  ...  M. 


\n 


(55) 


We  find  from  Eq.  (55)  that  the  RMS  phase  error  is  minimum  at  the  center  of  the  array, 

<rmin  «  0.234  A  tan  A<£  Ml/2,  degrees, 
and  maximum  at  the  edges  of  the  array 

‘’’max  —  2(1  —  0.75A/_I)(rmin. 

We  find  from  Eqs.  (56)  and  (57)  that  when 

(A,  M,  A<£)  -  (1%,  9,  45°), 
the  minimum  and  maximum  RMS  phase  errors  are 

O- min  -  0-5  -  0.70°, 


(56) 

(57) 


and 


O’ max  “  O’!  “  O-,  -  1.29°. 

Thus,  a  1%  RMS  error  in  the  power  splits  gives  rise  to  about  a  1°  RMS  phase  gradient  error. 


MODULES  WITH  SELF-COUPLING 


Our  analysis  thus  far  assumes  that  the  coupled  module  lattice  comprises  of  amplifier/limiter 
modules  each  configured  as  shown  in  Fig.  3.  However,  alternative  module  structures  are  conceivable, 
such  as  that  depicted  in  Fig.  12.  Comparison  of  Fig.  12  with  Fig.  3  shows  that  the  only  difference 
between  the  two  structures  is  the  presence  of  a  self-coupling  path  in  Fig.  12. 


14 


NRL  REPORT  8492 


A  transient  analysis  of  CAM  lattices  composed  of  self-coupled  modules  is  given  in  Appendix  F. 
It  turns  out  that  the  self-coupled  module.  Fig.  12,  has  less  desirable  performance  characteristics  than 
the  simpler  module  analyzed  previously. 

In  brief,  we  find  in  Appendix  F  that  (cf.  Eq.  (F34)) 

lA0„(r)  -  AM~>WT/r5  (58) 

for  self-coupled  modules.  Fig.  12.  The  settling  time  T,  in  Eq.  (58)  is  given  by 

r,  -  Tj j  1  +  y  r'  sec  Ad>|,  (59) 

where  r,  —  M2/5  is  the  settling  time  for  modules  without  self-coupling,  Eq.  (34).  Stable  equilibrium 
phase  distributions  exist  only  for  0  <  &<t>  <  w/2  when  F  >  0;  stable  distributions  exist  only  for  7r/2  < 
A<£  <  w  when  F  <  0.  It  follows  from  Eq.  (59)  that 

T,  >  r„  (60) 

always,  i.e.,  modules  without  self-coupling  settle  faster  than  modules  with  self-coupling. 

Figure  13  depicts  the  transient  development  of  A0|(t)  obtained  by  means  of  the  theory  developed 
in  Appendix  F,  for  the  same  parameter  values  as  Figs.  7  and  8.  Comparison  between  Figs.  13  and  8 
shows  that  the  phase  distribution  for  self-coupled  modules  evolves  more  slowly  with  time  and  is  tem¬ 
porally  smoother  than  the  phase  distribution  of  modules  that  lack  self-coupling. 

According  to  Eq.  (59), 

•  settling  is  faster  for  strong  coupling  (large  F)  than  for  weak  coupling  (small  F).  However,  settling 
time  cannot  be  decreased  indefinitely  by  increasing  F:  once  Fcos  A<t>  »  I  further  increases  in  F 
provide  no  further  improvement.  We  find  that 
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Fig.  13  —  Transient  development  of  A8|(t)  for  a  self-coupled  CAM 
lattice.  Assumed  parameter  values  are  (A#,  A<£0,  oM( 0))  —  (9,  60°. 
45°),  for  all  curves.  Strong  coupling,  F  ~  10,  results  in  faster  settling 
than  unity  coupling,  F  “  I  Solid  curves  obtained  by  numerical  solu¬ 
tion  of  Eq.  (FI);  dashed  curves  obtained  by  one-mode  approxima¬ 
tion,  Eq.  (F34).  Solid  and  dashed  curves  for  F  “  10  are  too  close  to 
distinguish  on  this  plot. 


lim  T,  -  r,.  (61) 

f— «• 

i.e.,  we  can  do  no  belter  ihan  the  seilling  lime  obtainable  by  eliminating  the  self-coupling.  Also, 

•  settling  time  becomes  progressively  longer,  approaching  infinity,  as  A<f>  approaches  90°. 

The  longer  settling  time,  Eq.  (60),  and  the  dependence  of  settling  time  on  phasing.  Eq.  (59),  are 
undesirable  attributes  that  lead  us  to  favor  modules  without  self-coupling  over  modules  with  self¬ 
coupling. 

TWO-DIMENSIONAL  ARRAYS 

We  now  briefly  discuss  how  our  results  for  1-D  CAM  lattices  can  be  extended  to  2-D  structures, 
appropriate  as  feeds  for  2-D  phased  arrays. 

We  adopt  an  abbreviated  convention  for  depicting  CAM  lattices,  as  illustrated  in  Fig.  14  for  the 
previously  analyzed  1-D  structure.  The  empty  circles  in  Fig.  14  denote  amplifier/limiter  modules  ter¬ 
minated  by  antennas;  filled  circles  denote  modules  with  nonradiating  terminations.  Thus,  Fig.  14  is 
equivalent  to  Fig.  2. 


•xp  JO  »xp  I  sa* 

Fig.  14  —  One-dimensional  CAM  lattice.  Empty  circles  denote 
amplifier/limiter  modules  terminated  by  an  antenna;  filled  circles  denote 
modules  with  nonradiating  terminations.  This  figure  is  equivalent  to  Fig 
2 
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f  igure  15  is  a  2-L)  generalization  of  Fig.  14.  Based  on  results  already  given,  it  is  easily  shown  that 
the  equilibrium  phase  distribution  of  the  Fig.  15  lattice  is  the  desired  2-D  linear  phase  gradient.  We 
can  see  by  inspection  that  the  settling  time  for  the  2-0  square  lattice  must  satisfy  the  constraint 

ts(I-0)  <  rs(2 -D)  <  2t,( I -D). 

where 

ts  (1-D)  —  A/2/5  is  the  settling  time  of  1-D  edge  lattice, 

and 

rs  (2-D)  is  the  settling  time  of  2-D  lattice. 


Fig.  15  —  Two-dimensional  CAM  lattice  suitable  for  excit¬ 
ing  a  2-D  phased-array.  While  only  three  phase-shifters  are 
required  to  excite  a  rectangular  structure  of  arbitrary  size, 
settling  time  increases  as  the  square  of  the  largest  linear 
dimension  of  the  lattice.  Thus,  additional  phase-shifters 
may  be  required  to  satisfy  bandwidth  requirements. 
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Appendix  A 

THE  PHASE  DIFFUSION  EQUATION 


Taking  the  imaginary  part  of  Eq.  (4)  we  obtain 

sin  9„(t  +  1)  -  /Vb(t)  ■  S„(t).  (A I > 

where,  by  definition, 

S„(t)  -  Bn  sin  0„+\(t)  4-  Fn  sin  (A2) 

Equation  (A2)  may  be  rearranged  to  obtain, 

S„(t)  ”  y(#„  +  /^Ksin  0„+ 1  +  sin  0„_,)  +  -j(£„  -  F„)(sin  9„+[  -  sin  0„_,).  (A3) 

Applying  some  standard  trigonometric  identities  to  Eq.  (A3), 

S„(t)  -  cos  y(0„+i  +0*_i)  cos  -«,-)) 

i(Bn  +  Fn)  tan  y(0fl+l  +  +  ( Bn  -  F„)  tan  |(@„+I  -  «„_,)}.  <A4> 


Similarly,  taking  the  real  part  of  Eq.  (4)  we  obtain 

cos  9„(t  +  1)  -  (V»(t)  •  C„(r). 


where,  by  definition, 

C„(t)  -  Bn  cos  (r)  +  F„  cos 

Equation  (A6)  can  be  recast  as 

C„(t)  -  cos  y  (flB+i  +  »„-[)  cos  y(0„+i  -  0„-i) 

{{Bn  +  F„)  -  ( Bh  -  F„)  tan  y(0„+,  +  0n-\)  tan  y  (9n+l 


(A5) 

(A6) 


0„-,)).  (A7) 


Dividing  Eq.  (Al)  by  Eq.  (A5),  and  substituting  Eqs.  (A4)  and  (A7)  for  S„(t)  and  C„(t), 


tan  0„(r  4-  1) 


tan  y(0„+,  +  0„-|)  -  y  e„  tan  y(0(t+i  -  9„_ ,) 


1  +  y  «*  tan  y(0„+)  -  #„_,)  tan  y(0„+)  +  0*-|) 


where 


e-  =  2 


B„ 


Fn+  Bn 

We  see  immediately  from  Eq.  (A8)  that,  for  F„  =* 

I 


#fl(r  +  1)  -  y [»„+, (r )  +  0„_, (r)J,  «„  -  0. 


(A8) 


(A9) 


(A10) 
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Subtracting  0„(t)  from  each  side  of  Eq.  (AIO), 

[ttn(T  +  1)  -  0„(t)\  -  yfe»„  +  ,(T)  -  20„(t )  +  0„_,(r) J.  (AID 

Equation  (All)  has  the  form  of  a  one-dimensional  heat  flow  equation  in  which  the  lime  and  space 
derivatives  have  been  replaced  by  finite  differences.  We  refer  to  Eq.  (All)  as  the  phase  diffusion  equa¬ 
tion  for  perfect  modules,  i.e.,  modules  for  which  =  0. 


More  generally,  further  simplification  of  Eq.  (A8)  when  e„  ^  0  requires  that  we  assume 


y«.  ‘an  j(»„+l  «  I, 


so  that 


1 


fan  -  0„_|)  =  tanjy  e„  tan  y (»„+,  - 


I 


2  "  2 

It  follows  from  Eqs.  (A8),  (A  13),  and  the  trigonometric  equality 

tan  A  -  tan  B 


tan  (A  —  B) 


1  +  tan  A  tan  B  ’ 


that 


<Ur  +  1) 


(A  12) 


(A  13) 


#ff+i(r)  +  0B-| (t)  -  e„  tan  y[0„+l(r)  -  0„_,(r)J}. 


(A14) 


We  refer  to  Eq.  (A  14)  as  the  general  phase  diffusion  equation. 
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Appendix  B 

EQUILIBRIUM  PHASE  DISTRIBUTIONS 


ANALYSIS 

For  a  given  initial  phase  distribution,  Eq.  (20),  and  a  specified  step-increase  in  the  boundary  phas¬ 
ing  at  r  =  0 

aM(0)  =  mod  :,(«£- 0J,).  (Bl) 

the  CAM  network  may  evolve  either  of  two  possible  phase  distributions  (cf.  Fig.  5).  We  presently 
address  the  question  of  which  of  these  two  final  phase  distributions  is  ultimately  established. 

We  have  found  by  computer  simulation  that  the  parameter  s  in  our  analytically  derived  equili¬ 
brium  phase  distribution,  Eq.  (19),  can  be  determined  by  inspecting  the  phase  distribution  after  just  a 
single  delay  time,  r  =  l.  Only  0w_i,  one  module  removed  from  the  point  of  disrupted  excitation,  has 
changed  from  its  initial  value  at  this  time.  Defining  the  time  gradient  of  phasing  as 

a„(r)  =  mod2ir{e„(r)  —  0„(r  -  1)},  (B2) 

we  find  that 

=  m°d2W{fl*/_i(l)  —  #a#-i (0))  n  s  m  j.  (B3) 

The  meaning  of  Eq.  (B3)  is  interpreted  with  the  aid  of  Figs.  Bl  and  B2.  Assuming  that  the  phase 
distribution  is  given  by  Eq.  (20)  at  time  r  <  0,  we  must  determine  whether  0^,(1)  is  greater  or  less 
than0w_!(O).  First,  we  note  that 

0„(O)  -  0„(-l),  n  -0,  1,  2.  ...  (A/-  1).  (B4) 

and 

0„(1)  -  0,(0),  n  -  0,  1,  2,  ...  (A/-  2).  (B5) 

since  the  phase  disturbance  initiated  at  ( n ,  r)  —  (M,  0)  propagates  at  the  rate  of  one  module  per  unit 
of  normalized  time.  Comparing  Fig.  Bl  and  Fig.  5,  we  see  that  Fig.  Bl  corresponds  to  an  incipient 
s  —  0  phase  distribution.  Similarly,  comparison  with  Fig.  5  shows  that  Fig.  B2  corresponds  to  an  inci¬ 
pient  s  —  1  distribution. 

Setting  (n.  t)  =*  ( M  —  1,  0)  in  Eq.  (4)  we  obtain 

«,/«-!<»  »  Nw_,(0)Fw_ ,[<-'* w'0*  +  (B6) 

Multiplying  both  sides  of  Eq.  (B6)  by  exp  (— y0w_|(O)j  we  obtain 

.  Aiw_1(0)FM_,le'Mw<0’  +  e-;4’*-'10’],  (B7) 

where  aM_|(l)  is  defined  by  Eq.  (B2)  and  the  quantities  A0„(t)  are  defined  by 

A0„(t)  =  0„(r)  —  0„_,(t).  (B8) 

From  Eqs.  (B4)  and  (B8) 

A0M_,(O)  -  0W_,(-1)  -  0w_2<-D.  (B9) 

From  Eqs.  (B9)  and  (20), 

A0W_,(O)  -  (B10) 
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«n(1) 


Fig.  B1  —  Incipient  5-0  equilibrium  phase  distribution.  Circle  (O)  denotes 
value  it  t  ■  0.  Cross  (x)  denotes  value  at  r  "  1.  Triangle  (A)  denotes 
values  that  are  the  same  at  r  ~  0  and  r  —  1. 


Fig.  B2  —  Incipient  S  “  I  equilibrium  phase  distribution. 
Symbols  O,  *,  and  A  as  in  Fig.  BI. 

Similarly,  we  can  show  from  Eqs.  (5b),  (B4),  (B8),  and  (20)  that 


From  Eqs.  (Bl)  and  (5b) 


A0jy(O)  “  -F  M( A<^  —  A0q). 


(BID 


a*y(0)  “  mod2„(®)&  ®m)  “  A/(A$  —  A<6(JI.  (BI2) 

Substituting  Eqs.  (BIO)  and  (BID  into  Eq.  (B7),  and  making  use  also  of  Eq.  (B 1 2) ,  we  obtain 

eJ°M~ -  A'M_|(0)F’w_l[^4*V’M<m  +  e-'A4°].  (B13) 

Equation  (B13)  is  given  a  phasor  diagram  representation  in  Fig.  B3,  for  positive  coupling  (Fm~\  >  0). 


We  see  from  Fig.  B3  that  «w_|(l)  <  jt  so  long  as 

0  <  «M(0)  <  (tt  —  2A<£„).  (B14) 

In  interpreting  Eq.  (B14)  we  recall  that  aw(0)  is  the  change  in  edge  phasing  imposed  by  (he  phase- 
shifter,  as  given  by  Eq.  (BI2).  If  we  do  not  increase  the  edge  phasing  by  too  large  an  amount,  Eq. 
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Fig.  B3  —  Phasor  diagram  derived  from  Eq.  (B13),  assuming  F>  0.  The 
5  —  0  phase  distribution  (cf.  Fig.  5)  evolves  when  angle  aw_ ,0)  <  ir, 
i.e.,  for  0  <  aw(0)  <  (v  -  2A4>q) 


(BI4)  will  be  satisfied  and  the  s  —  0  phase  distribution  will  evolve  after  a  period  of  time  (cf.  Fig.  5). 
However,  if  the  edge  phasing  is  increased  to  such  an  extent  that  Eq.  (B14)  is  not  satisfied,  we  find 
instead  that  the  s  =  1  phase  distribution  will  ultimately  evolve. 


Figure  B3  and  the  conclusions  derived  from  it  pertain  to  positive  coupling  (F  >  0).  The 
corresponding  phasor  diagram  for  negative  coupling  (F  <  0),  also  derived  from  Eq.  (B13),  is  shown  in 
Fig.  B4. 

In  general,  the  equilibrium  phase  distribution  can  be  predicted  from  Eq.  (19)  and  the  following 
set  of  rules  for  assigning  a  value  to  the  parameter  s. 


For  positively  coupled  modules  (F>  0;  0  <  A<fo>  <  w/2), 

0,  0  <  aw(0)  <  (ir  —  2A<£q) 
1 ,  otherwise. 

For  negative  coupling  (F  <  0;  ir/2  <  A 0O  <  ir), 

0,  0  ^  oijk(0)  <  (3ir  —  2A0o) 
1 ,  otherwise. 


(B15) 


(B16) 


Fig.  B4  —  Phasor  diagram  derived  from  Eq.  (B13),  assuming  F  <  0. 
Note  that  >  90*.  The  5  —  0  phase  distribution  (cf.  Fig.  5) 
evolves  when  <*w_|(l)  <  ir,  i.e.,  for  0  <  au(0)  <  (3ir  -  2A4q)- 
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PHASING  UP  (s  =  0  Phase  Distributions) 

By  "phasing  up"  we  mean  the  proeess  of  modifying  the  edge  phasing  to  achieve  a  new  equilibrium 
phase  gradient  that  is  larger  than  the  initial  gradient,  AW„  >  A«4n-  Inspection  of  Pig.  5  shows  that  phas¬ 
ing  up  corresponds  to  the  establishment  of  an  s  =  0  phase  distribution. 

From  Eqs.  (19).  (21),  and  (23)  (assuming  F>  0), 

AO  „  =  A<£o  +  [a  m(0)/ M],  (BI7a) 

subject  to  the  requirement  that 

0  ^  aw(0)  <  (it  -  2A*0>.  (B17b) 

where  A <t>u  is  the  initial  phase  gradient,  A0„  is  the  new  phase  gradient,  and  a ,4/(0)  is  the  step-change  in 
edge  phasing.  According  to  Eq.  (B17),  the  new  phase  gradient  is  proportional  to  the  step-change  in 
edge  phasing,  aju(0 ),  so  long  as  a,w(0)  is  not  too  large.  While  Eq.  (B17)  is  applicable  only  to  positively 
coupled  lattices  ( F  >  0),  an  analogous  formulation  is  readily  derived  for  F  <  0.  The  upper  limit  on 
A0„  established  by  Eq.  (B17)  is 


A 0„  <  [A4>0  +  (w  —  2A<(>0)/M].  (BI8) 

For  example,  we  assume  that 

(A*0.  M)  -  (60°,  9),  (B19) 

i.e.,  that  we  have  an  8-module  feed  initially  phased  to  60°  per  module.  If  we  step-increase  the  edge 
phasing  by  just  less  than  the  upper  limit  indicated  in  Eq.  (B17b),  say, 

aM(0)  *=  0.99(tt  -  2A*0)  -  59.4°,  <B20> 

the  lattice  will  eventually  settle  into  a  new  phase  distribution 

A0„  =  60°+  (59.4 °/9)  “  66.6°,  (B21) 

just  less  than  the  upper  limit  established  by  Eq.  (B18),  A0„  <  66.67°.  If  we  attempt  to  increase  A 9„ 
above  66.67°  by  step-increasing  0M  by  an  amount  aM(0)  >  (ir  —  2A <£0),  the  phase  gradient  decreases 
rather  than  increases,  as  discussed  below. 

PHASING  DOWN  (s  *  1  Phase  Distributions) 

By  "phasing  down"  we  mean  the  process  of  modifying  the  edge  phase  to  achieve  a  new  equili¬ 
brium  phase  gradient  that  is  smaller  than  the  initial  gradient,  A 0n  <  A<*>„.  Phasing  down  corresponds  to 
the  evolution  of  an  s  =  I  phase  distribution  (cf.  Fig.  5). 

From  Eqs.  (19),  (21),  and  (23)  (assuming  again  that  F>  0), 

A0„  =  A<£q  —  127 r  -  aM(0))/M,  (B22a) 

subject  to  the  requirement  that 

(it  —  2A</>0)  <  ^  277,  (B22b) 

where  A<£„,  A 0„,  and  a,w(0)  arc  defined  as  in  Eq.  (BI7).  The  lower  limit  on  A0„  established  by 
Eq.  (B22)  is 

A0„  >  (Ad>o  -  ( 77  +  2A<£(|)/ M).  (B23) 
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For  example  we  assume  (A<6 0,  M)  =  (60°,  9),  as  in  Eq.  (B19),  and  lhai  the  edge  phasing  is 
step- increased  by 

a ,t/(0)  =  1,01  (it  -  2A <*>0)  =  60.6°  (B24) 

We  find  from  Eqs.  (B19),  (B22a),  and  (B24)  that  the  lattice  eventually  establishes  a  uniform  phase  gra¬ 
dient 

A0„  -  60°+  (60.6°-  360°)/9-  26.73°,  (B25) 

just  greater  than  the  lower  limit  obtained  from  Eq.  (B23),  A0„  >  26.67°. 

We  have  seen  that  step-increasing  the  edge  phase  9U  by  59.4°  increases  the  phase  gradient  from 
60°  per  module  to  66.6°  per  module  (cf.  Eqs.  (BI9),  (B20),  and  (B2D).  A  slightly  larger  edge  phase 
increase  of  60.6°  decreases  the  phase  gradient  from  60°  per  module  to  27°  per  module  (Eqs.  (B24)  and 
(B25)). 


PHASING  LIMITS  AND  MULTISTEP  PHASING  PROCESSES 


Starting  from  an  initial  phase  gradient  of  A <£0  per  module,  it  follows  from  Eqs.  (BI8)  and  (B23) 
that  there  is  no  step-change  in  edge  phase  that  will  develop  a  phase  gradient  outside  the  range 

(A0O  -  (w  +  2A<6o)/A/1  <  A 0n  <  [Ad>0  +  (n  -  2A^0)/Wl.  (B26) 

Assuming  for  example  that  (A 4>0,  M)  =  (60°,  9)  we  obtain 

26.67  °  <  A0„  <  66.67°.  (B27) 

This  apparent  limitation  on  the  range  of  A 0„  is  actually  readily  overcome.  A  multistep  process  can  be 
used  to  sweep  the  phase  gradient  outside  the  range  indicated  by  Eq.  (B26). 


Considering  the  phasing  up  process,  for  example,  we  generalize  Eq.  (B17)  as  follows: 

A e'ni+n  -  Ml*  +  lati+))/M),  (B28a) 

subject  to  the  requirement  that 

0  <  a#+l)  <  frr  -  2A »ln),  (B28b) 

where 

A0<°>  =  A*„,  (B28c) 

and  where  A#^'1,  /  =  1,  2,  . . .  ,  is  the  equilibrium  phase  gradient  that  evolves  after  the  Ah  step- 
increase  in  edge  phasing,  a|}’. 


Continuing  our  previous  example,  we  find  from  Eqs.  (B28),  (BI9),  (B20),  and  (B21)  that 

59.4°| 

A0‘"  -  66.6°j- 


(B29) 


A  new  phase  gradient  A0*2’  >  A0*n  may  be  established  by  initiating  another  step-increase  in  edge 
phasing,  so  long  as  we  observe  the  constraint  imposed  by  Eq.  (B28b).  Choosing,  for  example, 


«J/2)  -  0.99[rr  -  2A0(<I)]  -  46.33° 


(B30) 


we  establish  a  new  equilibrium  phase  gradient 

A0^2)  -  66.6 6  4-  (46.33  79)  -  71.75°, 


(B31) 


once  the  transient  has  decayed.  The  process  may  be  continued,  and  ever-increasing  phase  gradients 
thereby  achieved. 
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We  obtain  the  multisiep  prescription  for  phasing  down  by  generalizing  Eq.  (B22), 

-  [2jt  -  a^"]/M,  (B32a) 

subject  to  the  requirement  that 

l»r  -  2Afl‘'M  <  «i}+"  <  2tt.  (B32b) 

where  A0J'1  and  a^+u  are  defined  as  in  Eq.  (B28). 

The  maximum  increase  in  the  phase  gradient  when  phasing  up  is  computed  from  Eq.  (B28)  as 

[A0i/+u  -  Afl^'l  <  [it  -  2A0',>1/A/.  (B33) 

The  maximum  decrease  in  the  phase  gradient  when  phasing  down  is  obtained  from  Eq.  (B32)  as 

(A»^‘  -  A6><,+1,1  <  In  +  2A»‘/,)/M  (B34) 

Comparing  Eqs.  (B33)  and  (B34)  we  see  that  phasing  down  the  lattice  generally  requires  fewer 
steps  than  phasing  up.  Considering  our  previous  example,  Eq.  (B19),  we  see  that  the  lattice  can  be 
phased  down  (60°  -  26.7°)  =  33.3°  in  a  single  step.  By  contrast,  a  single  step-change  in  the  edge  phas¬ 
ing  can  phase  up  the  gradient  by  no  more  than  (66.6°  -  60°)  =  6.6°. 

An  important  conclusion  derivable  directly  from  Eq.  (B33)  is  that  as  AO^'  approaches  90°  per 
module,  lA0,J,+,>  —  Afl^1]  approaches  zero,  and  no  further  increases  in  the  phase  gradient  are 
achievable.  Thus,  phase  gradients  90°  <  Afl„  <  180°  are  unattainable  for  positively  coupled  lattices. 

Our  discussion  following  from  Eq.  (B17)  pertains  only  to  positively  coupled  lattices,  F  >  0.  A 
parallel  formulation  for  F  <  0  shows  that 

•  phase  gradients  0°  <  \0„  <  90°  are  unattainable  with  negatively  coupled  lattices,  and 

•  negatively  coupled  lattices  can  be  phased  up  over  the  range  90°  ^  A 0„  <  189°  much  more  rapidly 

than  they  can  be  phased  down. 
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Appendix  C 

TRANSIENT  ANALYSIS 


In  this  appendix  we  derive  an  exact  solution  to  the  equation 

9„(t  +  l)  -  y[fl,+1(T)  +  »„_i(t)1,  (Cl) 

subject  to  the  initial  conditions 

0„(t)  -  wAdo*  h  ■*  I.  2.  ...  M.  t  <  0,  (C2) 

and  boundary  conditions 

0o(t)  -  0 

0"(t)  -  MA9„(°») 

Also,  the  equilibrium  solution  of  Eq.  (Cl)  in  the  limit  of  infinite  time  is  given  by 

0„(«>)  -  /jA0„(°o),  (C4) 

where  A0„(<»)  is  a  constant  independent  of  n. 

We  begin  our  solution  of  Eq.  (Cl)  by  transforming  to  a  new  phase  variable  <|»„(t),  subtracting  the 
known  equilibrium  value  from  9„(r ), 

i Mr)  =  0n(r)  -  9„(oo)  -  0n(r)  -  »A0„(°°).  (C5) 

Thus,  <9„(t)  represents  just  the  transient  part  of  0„(r), 

lim  i J>b(t)  -  0-  (C6) 

r— oo 

Applying  Eq.  (C5)  to  Eq.  (Cl)  we  obtain  an  identical  difference  equation  for  0„(r), 

0„(t  +  1)  -  -y-(d*„+i (t)  +  (r)J.  (C7) 

The  new  boundary  and  initial  conditions  are  obtained  from  Eqs.  (C2),  (C3),  and  (C5)  as 

</»0(t)  -  ^m(t)  —  0,  (C8) 

and 

«„0  =  I K(0)  -  n[\<t>Q  -  A0„(~)J,  n  -  1,  2  ,  ..AM  -  I).  (C9) 

The  point  of  our  variable  transformation,  Eq.  (C5),  is  to  provide  us  with  homogeneous  boundary  con¬ 
ditions,  Eq.  (C8),  after  which  the  initial  value  problem  becomes  solvable  by  separation  of  variables. 

We  look  for  solutions  to  Eq.  (C7)  of  the  form 

<Mt)  “  Z„/(t).  (CIO) 

Substituting  Eq.  (CIO)  into  Eq.  (C7)  we  obtain 

z„/(T  +  n-  j(z„+1  + z„_,)/(t). 


Thus, 


/(t  +  1  )//(r)  -  (Z„+l  +  Z„-\)/2Z„  -  X. 


(CIO 
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where  the  parameter  A  is  dependent  on  neither  r  nor  n.  From  Eq.  (Cl  1) 

/(r  +  I)  =  X/(t).  (Cl  2) 

Equation  (Cl 2)  is  solved  either  by  inspection  or  by  standard  difference  equation  methods  to  obtain 

fir)  =  /(0)Xr,  (CI3) 

where  /' ( 0 )  is  an  undetermined  constant. 


Equation  (Cl  1)  also  provides  us  with  an  algebraic  eigenvalue  equation  for  Z„, 
j(Zn+1  +  Z„_,)  =  XZ„,  »-l.  2.  ..AM  -  1), 
where,  from  Eqs.  (C8)  and  (CIO), 


(Cl  4) 


(C15) 


20  =  ZM  -  0. 

As  may  be  verified  by  substitution,  the  eigenvalues  \k  and  corresponding  eigenvectors  Z„u>  of  Eq. 
(C14)  are 


X.k  —  cos  ( kir/M ), 


and 


Zn(*'  =  (2/ A/) 1/2  sin  (wkn/M), 
for  k  =  I,  2,  ..AM  -  1).  Eigenvectors  (Cl  7)  satisfy  the  orthogonality  condition 


Y  Z}k'Z“'  =  8U  ■ 
1 


1,  k  -  / 
0 ,  k  *  l. 


The  general  solution  of  Eq.  (C7)  is  obtained  from  Eqs.  (CIO)  and  (Cl 3)  as 

<Mr)  -  Y  A*  Z"U>  x*’ 


(Cl  6) 
(C17) 

(Cl  8) 

(C19) 


where  the  constants  Ak  must  be  determined  from  the  initial  conditions.  We  set  r  =  0  in  Eq.  (C19), 
multiply  by  Zff<M,  and  sum  over  index  n  to  obtain 

f«-t 


M- 1  M- 1 

I  *,(0)z:n  -  £ 

n—  I  A“  f 


I  zrz:1 

n—  I 


r<*> 


Thus,  from  Eqs.  (Cl 8)  and  (C20), 

M- 1 

-4*  "  X  ^ »° 

n*  1 

Substituting  Eq.  (C9)  for  i|»„0  and  Eq.  (Cl 7)  for  Z„ul  into  Eq.  (C21)  we  obtain 

M-l 

Ak  -  (2/A/)1/2[A0„(°o)  -  £  «  sin  (nkn/M). 

n—\ 

The  summation  in  Eq.  (C22)  can  be  performed  analytically,  with  the  result 

Ak  -  ( M/2)'/2lM„(°o )  -  A«0J(-l)*cot(7rX/2M). 
From  Eqs.  (C19),  (C23),  (C16),  and  (C17), 


where 


<t)  -  iV:*’(t). 

A- 1 


■u,(r)  -  lA0„(«>)  -  A<^0H-l)*cot|-^g 
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sinl 


nkn 


M 


cos 


X7T 


(C20) 

(C21 ) 

(C22) 

(C23) 

(C24) 

(C25) 
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Since 


—  cos 


(M  -  k)7T 
M 


we  see  from  Eq.  (C25)  that  i has  the  same  time  dependence  as  i This  suggests  writing 
Eq.  (C24)  in  the  form 


<I>„(t)  =  j  +  >PnM  *  *  (r)), 

k-1 


where  [AY/2 ]  is  the  integer  part  of  A//2, 


[M/2]  =  | 

We  can  show  from  Eqs.  (C25)  and  (C26)  that 


(M/2)  ,  M  even 
(M-D/2  ,  M  odd. 


•Mr)  =  2[A0„(oo)  -  Atfol  '  Jp'  i-lA&lS&m.  [cos  (kn/M) Y+*. 

fa  sin  (kn/M) 


where 


fl  =  y[I  +  (-l)n+r+w]. 


Equation  (C28)  is  the  exact  solution  of  Eq.  (Cl),  subject  to  Eqs.  (C2)  and  (C3),  the  initial  conditions 
and  boundary  conditions,  respectively. 


Appendix  D 

UNEQUAL  FORWARD  AND  BACKWARD  COUPLING 


FORMULATION 

The  equilibrium  phase  distribution  of  the  CAM  network  is  the  solution  of  Eq.  (13), 


where 


A0fl+1  -  A0„  +  e„  tan  y  (A +  A »„), 


_  1  ^  , 

e.  =  2  — - —  «  1. 


We  now  define  a  quantity  A<£„  that  measures  the  degree  to  which  the  actual  phase  gradient  A 9„  departs 
from  the  ideal,  i.e.. 


A*.  -  A 9n  -  A®“. 


From  Eqs.  (Dl)  and  (D3) 


where,  from  Eq.  (19), 


A*„+1  -  A*„  +  eB  tan  A0fl°  +  (A*„+l  +  A*B)  , 


A#?  -  AO  -  s(2r/M).  n  -  1,  2,  ...<♦/. 


If  the  CAM  network  is  to  be  useful  for  exciting  a  phased  array  antenna,  it  must  be  true  that 


y  (A***,  +  A*„)  «  A6>», 


from  which  we  can  show  that 


tan  A +  y  (A*„+l  +  A <f>„)  =  tan  A0°  +  y  (A*„+1  +  A*,)  sec2  A®J.  (D7) 


From  Eqs.  (D4)  and  (D7), 


A<^»+i  -  A tan  A 0°  +  y  «„(A*1I+,  +  A*„>  sec2  A0°. 


Since  development  of  an  exact  dosed  form  solution  to  Eq.  (D8)  appears  quite  difficult,  we  pro¬ 
pose  an  iterative  solution, 

A^V’  -  A*‘,+n  +  tan  A#n°  +  y  +  A<'>)  sec2  \0°„.  (D9) 

In  the  limit  of  fabricalionally  perfect  modules, 

lim  A0„  =  A»“.  (DIO) 

€  _  0 


The  ^ero-order  iterate  A0„m  is  thus  obtained  from  Eqs.  (D3)  and  (DIO) 

A<”  -  0. 


(Dll) 
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The  first-order  iterate  follows  from  Eqs.  (D9)  and  (Dll), 

+  «„  tan  A#®.  (D12) 

We  simplify  notation  by  dropping  the  superscript  on  A^1’  and  by  setting  s  -  0  in  Eq.  (D5)  to  obtain 

A4„-m  ”  A0„  +  c,  tan  A^.  (DI3) 


The  solution  of  Eq.  (D13)  is 


M 


A*,  -  K  -  tan  A*  £  n  -  1,  2,  ...  M,  (D14) 

as  may  be  verified  by  back  substitution  into  Eq.  (D13).  The  constant  K  in  Eq.  (D14)  is  independent  of 
the  module  index  n. 

Since  the  values  of  90  and  9M  are  imposed  by  the  phase  control  element  tcf.  Fig.  2),  they  are 
independent  of  whether  the  modules  are  fabricationally  perfect  or  imperfect;  thus, 


«o  -  «o° 

Defining 

boundary  conditions  on  tf>„  are  obtained  from  Eqs.  (DIS)  and  (D16), 

<t>o  "  4>m  *  0. 

We  impose  these  boundary  conditions  on  Eq.  (DI4)  to  determine  the  value  of  K. 

Toward  determining  the  proper  value  of  K  we  note  that 
M  M 

£  A^,  —  £  ($„  -  ^„_t)  -  (0a/  —  <£o)  “  0. 

/»— I  n- 1 

where  the  final  equality  in  Eq.  (D18)  follows  from  Eq.  (D17).  However,  from  Eq.  (D14), 

M  MM 

£  A*„  -  MK  -  tan  A*  £  £  «,. 

n-l  n—  I  p-n 


Thus,  from  Eqs.  (D18)  and  (D19) 


i  M  M 

tan  A <t>. 


(D  i  5) 

(D16) 

(DI7) 


(DI8) 


(D19) 


(D20) 


r—  1  p—r 


Substituting  Eq.  (D20)  into  Eq.  (D14)  we  obtain  the  general  first-order  approximate  solution  of 
Eq.  (D8), 


1  MM  Kt 

A<i"  ”  |77  II(r  £  «/>|  tan  A*,  n  -  1,  2.  ...  M. 

,  r-I  p-r  p-n 


(D21) 


We  obtain  the  solution  for  4>n  by  noting  that,  analogous  to  Eq.  (DI8), 

A <frr  “  —  0q)  “  4>n>  n  ™  )  •  2,  . . .  M. 


(D22) 
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From  Eqs.  (D2I)  and  (D22) 


**“  77  I  X‘p“  I  !«/-  tan  A*.  "  =  '•  2,  ...  M.  (D23) 

|  M  r-1  p-r  ,-|  p-,  J 

Equation  (D23)  correctly  reproduces  the  boundary  condition  at  n  =  M.  namely  d»  w  —  0.  However,  we 
should  not  be  worried  that  the  boundary  condition  at  n  —  0  is  not  reproduced  since  Eq.  (D22)  (from 
which  Eq.  (D23)  was  derived)  is  not  valid  for  n  -  0. 

Finally,  we  should  not  be  concerned  that  the  quantity  tM  appearing  in  some  of  the  above  equa¬ 
tions  is  not  strictly  defined,  since  terms  in  identically  cancel  from  our  final  results,  Eqs.  (D21)  and 
(D23).  Thus,  Eqs.  (D2I)  and  (D22)  can  each  be  written  in  a  form  in  which  «  M  does  not  appear: 


** 


i  M—  1  M- 1  M- 1 

X  X  *p~  X  V‘anA*  .  ,  n 

1  .  ,  n  -  M 

*- ,a" 4* 


M- 1  JM-I 


~\  X  X«P-I  XM‘anA* 


1  r 


r- 1 


.  n  -  1,  2.  ...  (M  -  1) 
,  n  -  M. 


The  latter  two  equations  evaluate  identically  to  Eqs.  (D21)  and  (D23),  respectively.  We  generally 
prefer  to  work  with  Eqs.  (D2I)  and  (D23)  for  obvious  reasons. 

CONSTANT  COUPLING  RATIO 

We  presently  assume  that 

«„  -€,  <i-l.  2,  ...  (M-  1),  (D24) 

i.e.,  we  assume  that  all  modules  in  the  CAM  network  are  identical,  even  though  they  are  imperfect. 
(We  recall  that,  by  definition,  a  perfect  module  is  one  for  which  F„  =■*  B„,  i.e.,  e„  =  0,  n  — 
l.  2.  ..AM-  1).) 


It  follows  from  Eq.  (D24)  that  the  summations  appearing  in  Eq.  (D2I)  are  readily  evaluated. 

HI--M 


jr  tp  -  (m  —  »  +  D«. 


Substituting  Eq.  (D25)  into  Eq.  (D21), 


where 


Atf*  =  -  ea„.  n  -  I,  2,  ...  M, 


a.  =  [n  —  (M  +  I  )/2|  tan  k<f>. 


(D26) 


(D27) 


We  see  from  Eq.  (D26)  that  the  error  in  the  intermodule  phasing  is  minimum  at  the  array  center  and 
maximum  at  its  edges. 

In  Fig.  II  we  compare  values  predicted  by  Eq.  (D26)  with  a  computer  solution  of  Eq.  (4),  for  the 
parameter  values 

(«,  A*.  M)  -  (10-2.  65°,  9).  (D28) 
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We  recall  that  Eq.  (D26)  is  a  first-order  approximation.  Later  in  this  appendix  we  calculate  a 
second-order  approximation  that  greatly  reduces  the  discrepancy  between  the  numerical  and  analytical 
calculations  evidenced  in  Fig.  1 1 . 

TAPERED  COUPLING  RATIO 

Perhaps  the  simplest  type  of  nonuniform  network  is  the  linear  taper,  for  which 

ne1(  n  -  1.  2,  ...  (M-  1).  (D29) 

From  Eqs.  (D29)  and  (D13), 

A^„+/  “  A0„  +  «t|  tan  Ad>.  (D30) 


It  can  be  shown  that  the  summations  appearing  in  Eq.  (D21)  evaluate  as  follows. 


M 


and 


[M{M  + 1)-  »<"-  »i. 

if  I  Ui/6)  (M+  I)  (2M  +  1). 


r—  1  p—r 

From  Eqs.  (Dll),  (D31),  and  (D32), 


\<t>„  -  n  -  1,  2,  ...  M. 


where 


b„  »  |^-J  I3«(ff  -  1)  —  ( M 2  -  1)]  tan  A<£. 

To  obtain  A0„  in  degrees  Eq.  (D33)  (like  Eq.  (D26))  must  be  multiplied  by  (180/ir). 

We  note  that  the  value  of  n  for  which  Eq.  (D33)  is  equal  to  zero  may  be  approximated  as 

«0  =  0.5  +  0.577  M,  A <f>„0  «-  0. 


(D31) 


(D32) 


(D33) 


(D34) 


(D35) 


In  Fig.  D1  we  compare  values  predicted  by  Eq.  (D33)  with  a  computer  solution  of  Eq.  (4),  for  the 
parameter  values 

(c„  A4,  M)  -  (0.0025,  65°,  9). 

The  predictions  of  the  approximate  analytic  solution,  Eq.  (D33),  appear  quite  close  to  the  exact  numer¬ 
ical  solution. 


CONSTANT  COUPLING  RATIO:  SECOND-ORDER  APPROXIMATION 


From  Eqs.  (D9)  and  (D24)  we  obtain 


Atfi+i  -  A0^J)  +  c  tan  A«°  +  j  eCA^V’i  +  A^'n)  sec2  A0°, 
where  A^J+|  is  the  first-order  approximation  given  by  Eq.  (D26).  From  Eqs.  (D26)  and  (D27) 


y(A0”,1  +A<^")-£ 


n  — 


M 


tan  A<£. 


(D36) 

(D37) 
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Fig  D1  —  Equilibrium  phase  error  for  an  8-module  feed  with  a  linearly 
tapered  imperfection.  Parameter  values  assumed  are  («|.  M )  — 
(0  0025,  65",  9).  Solid  curve  obtained  by  numerical  solution  of  Eq.  (4); 
dashed  curve  obtained  by  first-order  approximation.  Eq.  (D33).  Quantity 
is  defined  only  at  integral  values  of  n. 


Substituting  Eq.  (D37)  into  Eq.  (D36), 

A0„+|  ■*  A^„  +  c  tan  A^  +  e2|/i  -  -^J  sec2  At i>  tan  A <f>,  (D38) 

where  we  simplify  notation  by  dropping  the  superscript  from  A^2)  and  by  substituting  A <j>  for  A 
(i.e.,  setting  s  *■  0  in  Eq.  (D5)).  Equation  (D38)  may  be  written  as 

A^,+|  —  A +  E  tan  A4  +■  nE\  tan  A <t>.  (D39) 

where,  by  definition, 

E  =  * 1 1  —  y  t\ i  sec2  A<A 

£i  =  *2  sec2  A<£.  (D40) 

The  solution  of  Eq.  (D39)  may  be  written  by  superposing  the  solutions  of  Eqs.  (D13)  and  (D30). 
Thus,  from  Eqs.  (D26)  and  (D33), 

A*„  =  A*<2'  -  Ea„  +  (D41) 

where  a„  and  b„  are  given  by  Eqs.  (D27)  and  (D34),  respectively.  We  can  show  from  Eqs.  (D41)  and 
(D26)  that 

[A*<2>  -  A^n]  -  y{ln  -  (A/  +  l)/2]2  -  (Af2  -  l)/12)c2  sec2  A*  tan  A <f>.  (D42) 

As  shown  in  Fig.  11,  the  second-order  approximation  brings  our  analytical  results  much  closer  to  the 
exact  numerical  solution  of  Eq.  (4). 

RANDOM  COUPLING 

Rather  than  lake  e„  as  a  known  function  of  n  (c.g.,  Eqs.  (1)24)  and  (1)29)),  we  now  assume  that 
the  coupling  coefficients  F„  and  B„  are  independent,  identically  distributed,  random  variables.  The 
quantity  e„  is  thus  also  a  random  variable  as  a  consequence  of  Eq.  (1)2). 
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We  denote  the  mean  value  and  variance  of  a  random  variable  X  as  nx  and  cr2,  respectively.  Thus, 
by  assumption,  hf  =  Ms  and  a-t  =  trB,  while  nf  and  o-,  are  independent  of  the  module  index  n, 

fi(  -  £(ej,  n  -  1.  2.  ...  KM  -  I) 

<r,2- £{(«,, -M«>3).  n  -  I.  2,  ...  (A/-  1).  (D43) 

The  statistical  expectation  operator  is  denoted  as  £(•]  in  Eq.  (D43). 

Since  a  simple  relationship  exists  between  £„,  and  B„ ,  Eq.  (D2),  we  might  expect  simple  rela¬ 
tionships  to  exist  between  their  mean  values  and  variances.  These  relationships  are  derived  as  Eqs. 
(EU)  to  (E15). 

To  simplify  our  notation  we  define 

H„  =  £{A0„), 

P2  S  £{A*2},  (D44) 

cr*  =  {(A<£„  -  “  (P?  ~  M3)- 

An  expression  for  /*„,  the  mean  phase  error,  is  obtained  by  applying  the  expectation  operator  to  Eq. 
(D21).  Since  the  expectation  operator  commutes  with  the  summation  operations  we  find 

I  X><  -  i  M«j  tan  A*.  (D45) 

Since  the  random  variables  «„  are  identically  distributed,  Iheir  mean  values  /x,  are  all  identical  (Eq. 
(D43)).  Thus,  Eq.  (D45)  evaluates  precisely  as  for  the  constant  coupling  ratio  problem,  Eq.  (D24). 
As  a  simple  adaptation  of  Eq.  (D26)  we  find  that 

H„  -  n,\n  —  (M  +  l)/2]  tan  A<£,  n  -  1,  2.  . . .  (A/  -  I).  (D46) 

From  Eqs.  (D46)  and  (E14) 

H„  ■■  10-4v42{w  —  ( M  +  l)/2]  tan  A <t>,  (D47) 

where  A  is  the  percentage  fabricational  accuracy,  as  defined  by  Eq.  (El 3). 

For  example,  assuming  that 

(A,  M,  A<£)  -  (1%,  9,  45”),  (D48) 

we  find  the  maximum  expected  phase  error  is  just 

max  n„  *»  (Im  —  4  x  10-4  (rad)  —  0.023°.  (D49) 

n 

Thus  the  mean  phase  increment  is  very  nearly  a  constant, 

£{A0b}  -  fin  +  A<J*  =  A<*.  (D50) 

From  Eq.  (D50) 

E[9„)  =  nA<t>.  (D51) 

It  is  of  interest  to  determine  the  extent  to  which  we  may  expect  the  phase  distribution  A 9„  to 
deviate  from  its  mean,  Eq.  (D50).  For  this  purpose  we  need  to  evaluate  the  second-order  statistic 
£[A02]  -  £lA0„l2.  First,  we  note  that 

£{A02)  -  £(A02}  -  P2.  (D52) 

as  a  consequence  of  the  fact  that 

A 9„  -  A tf>„  4-  A 4>.  (D53) 
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We  recall  that  A#  in  Eq.  (D53)  is  deterministic  and  is  determined  by  the  boundary  condition,  Eq.  (5); 
i.e.,  A is  determined  by  the  phase  control  element  in  Fig.  2. 


From  Eq.  (D2I), 


P „  =  tan2  A6  •  P„, 


(D54) 


where 


Pn 


i  M  \l 

ill', 

m  r-l  p-r 


M 

I', 


Hi  M  M  M 

Ji  I  I**'  I« 


(D55) 


Four  terms  are  obtained  when  the  multiplication  inside  the  expectation  operator  in  Eq.  (D55)  is  per¬ 
formed.  Commuting  the  expectation  operator  with  the  summations,  we  find  from  Eq.  (D55)  that 


Pn  -  L  +  R„  -  Mn  -  £„, 

(D56) 

where 

■  M  M  M  M 

M  r-l  s-l  p-r  fc-s 

,},  (DS7) 

Rn  =  I  2 

p-n 

(D58) 

i  M  t  M  M 

".  =  i  I  1 1  £i-M 

m  s-1  ^-n  *-s 

(D59) 

and 

1  M  [  M  M  | 

m  r-l  1  p-r  k-n  J 

(D60) 

Following  our  assumption  that  ( F„ ,  B„)  are  statistically  independent  of  (£m,  B„ ),  for  n  ^  m,  we  find 
that 

fo-,2,  p  -  k 

£[«,€*!  -  <r,28„.*  -  |0  p^k  (D61) 

With  the  definition 

M  M 

S„V  =  <7,2  £  f>p,k. 

p~u  *-v 

(D62) 

it  follows  that  Eqs.  (D57)  to  (D60)  may  be  written  as 

i  M  M 

l-vS 

(D63) 

“  Snn, 

(D64) 

'■-ij s- 

(D65) 

and 

I  w 

E"  ”  A/  r5  Srn' 

(D66) 

It  follows  from  Eq.  (D62)  that  Sin.  =  £V1(.  Thus,  from  Eqs.  (D65)  and  (1366),  wc  see  that  M„  —  £„. 
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It  can  be  shown  from  Eq.  (D62)  that 

S„v  =  <t2[(M  +  1)  —  max(z/,  v)], 
where  niax(w,  v)  denotes  the  larger  of  uand  v. 


max(t/,  v) 

It  can  be  shown  from  Eqs.  (D63)  to  (D67)  that 


u  ,  u  >  v 
v  ,  v  >  u. 


I  -  o-,2  •  -7  M~'(M  +  l)(2 M  +  1), 
6 

R„  -  <t2(M  +  1  -  n). 


M„  -  Eh  -  <7 g  •  y  AT'lAHM  +  1)  -  n(n  -  1)]. 

Collecting  terms,  it  follows  from  Eqs.  (D54),  (D56),  (D69),  (D70),  and  (D71),  that 

PH  -  £[( A*,)2)  -  (<r2  tan2  A*)(l  +  Ar')|y  +  y|  -  n  +  n2(M  +  1)-'J. 


We  complete  the  square  on  n  in  Eq.  (D72),  to  obtain 


I  i'/2 

P„  =  o-,  tan  A <t>  -jy  M(\  -  AT2)  +  XT'  In  -  (A/  +  1)/2]2|  . 
Equation  (D73)  clearly  displays  the  expected  symmetry  about  n  =  (M  4-  l)/2. 

From  Eqs.  (D73)  and  (E15),  and  multiplying  by  <  1  80/-jt ) ,  we  obtain 

P„  -  -7—  A  tan  A^l-i-  A/Cl  -  AT2)  +  M~'ln  -  (M  +  1)/2]2|  ,  degrees. 


(D67) 

(D68) 

(D69) 

(D70) 

(D71) 

(D72) 

(D73) 

(D74) 


Multiplying  Eq.  (D47)  by  (180/7r), 

fi„  -  -jyr  A2  tan  A <filn  -  (M  +  l)/2),  degrees.  (D75) 

Froms  Eqs.  (D44),  (D74),  and  (D75) 

*1  -  (Pi  -  pD  -  Pi.  (D76) 

so  long  as 

M  A2  «  2  x  IO4.  (D77) 

Assuming  that  Eq.  (D77)  is  valid,  and  also  that 

AT2  «  1,  (D78) 

we  find  from  Eqs.  (D74),  (D76),  (D77),  and  (D78)  that 

(7„  =  {£[(A0„)21  -  [£(A0„)I2}l/2  (D79) 

<t„  -  0.234  A  tan  A<Z>A/I/2{I  +  UM~2[n  -  (M  +  l)/2]2)'/2,  degrees,  n  -  1.  2,  ...  M.  (D80) 

We  recall  that  A  in  Eq.  (D80)  is  the  percentage  fabricational  accuracy,  defined  by  Eq.  (El 3);  (A/  -  I) 
is  the  number  of  modules;  A/A<£  is  the  phase  shift  established  by  the  phase  control  element  in  Fig.  2;  n 
is  the  module  index;  <r„  is  the  root-mean-square  (RMS)  deviation  of  A0„  —  9„  -  0„_|  from  its  desired 
value,  A0„  -  A<£. 
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Appendix  E 

RELATIONSHIPS  BETWEEN  COUPLING  STATISTICS 
AND  POWER  SPLITTER  ACCURACY 


Our  objective  in  this  appendix  is  to  derive  expressions  for  the  mean  and  variance  of 


e 


n 


2 


F„~Bn 

Fn  +  Bn 


(El) 


in  terms  of  the  mean  and  variance  of  the  power  split  coefficients  F„  and  Bn. 

As  a  matter  of  notation,  the  expectation  operator  is  denoted  as  £[■].  The  mean  and  variance  of 
random  variable  X  are  denoted  as  fix  and  cr 2,  respectively.  We  assume  that  random  variables  F„  and 
B„  are  independent  and  identically  distributed.  Thus 


Mr  ”  F-b 
<t  f=*  cr  a, 

for  all  values  of  n. 


(E2) 


Our  results  are  obtained  as  a  special  case  of  a  theorem  appearing  on  p.  14]  of  Ref.  El.  We  state 
the  general  form  of  the  theorem  first. 

If  the  random  variable  Z  is  a  function  of  the  two  random  variables  £and  fl, 

Z  =  H(F,  B),  (E 3) 

the  mean  and  variance  of  Z  may  be  approximated  by  the  following: 

Fz  —  H(nF,  Fb)  +  y  [(bFFH)  a-}  +  (9 is//)  <r£]  (E4) 

07  =  (dFH)2  <rf  +  (bBH)'  <r|.  (E5) 

a 

The  partial  derivatives  9f  -  — ,  etc.,  in  Eqs.  (E4)  and  (E5)  are  evaluated  at  ( F B)  —  (p.r,  nB). 

or 


For  the  special  case 

Z  -  F  B~\ 

and  making  use  of  Eq.  (E2),  we  find  from  Eqs.  (E4)  and  (E5)  that 

F’z  —  1  +  (o-jt/mf)2. 

and 

o-j2  =  2  ((rf/Mf)2- 

However,  from  Eqs.  (El)  and  (E6) 

F,-Fz~  1. 
and 

„  2  _  „  2 

(J  t  (T.. 


(E6) 

(E7) 

(E8) 

(E9) 

(E10) 
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From  Eqs.  (E7)  through  (E10). 


—  itTf/pf)2. 

(Ell) 

and 

(crF/pF). 

(E12) 

If  we  define  the  fabricational  accuracy  A  as 

A  -  100  {<r f/p  f)  (percent). 

(E13) 

Eqs.  (Ell)  and  (E12)  become 

p,  -  10~4  A2, 

(E14) 

and 

a,  -  y/i  X  1  Qr2A. 

(E15) 

For  fabricational  accuracies  A  = 
magnitude  larger  than  p,.  Thus, 
variable. 

1%,  we  see  from  Eqs.  (E14)  and  (El  5)  that 
we  are  often  able  to  ignore  p„  regarding  i 

<r,  is  about  two  orders  of 
e  as  a  zero-mean  random 

jcr  ,J  .  P  “  k 

Ek,**]  -  crXk  “  Jo  f  p  *  k 

(D61) 
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Appendix  F 

MODULES  WITH  SELF-COUPLING 


We  can  show  that  the  transient  phase  distribution  of  modules  with  self-coupling,  Fig.  12,  evolves 
according  to  the  equation 

-  /V„(t) [</»"’  +  Be'*"* |,r’  +  FeS”-'"),  (FI) 

analogous  to  Eq.  (4).  We  simplify  the  analysis  in  the  remainder  of  this  appendix  by  assuming  equal 
forward-  and  backward-coupling. 


With  the  definitions 


F  -  B. 

A0„(t)  —  0„(t)  —  0„_|(t), 

1 V  =  Nn  (t  ) , 

a  =  9„(t  +  1)  —  0„(t), 

A  =  y[Aflff+,(r)  +  A<Ut)], 
S  =  ylA0„+l(r)  -  A0„(r)l, 


Eq.  (FI)  may  be  written  as 


e>a  -  Af(l  +  IFe?  cos  A). 

Taking  the  real  and  imaginary  parts  of  Eq.  (F8)  we  obtain 

sin  a  -  2 NF  cos  A  sin  S, 


cos  a  -  N(\  +  IF  cos  A  cos  S). 


Dividing  Eq.  (F9)  by  Eq.  (F10), 


2  F  cos  A  P 

tan  a  -  -  —  — - - - r  sin  S. 

1  +  2F  cos  A  cos  S 


Inspection  of  numerous  numerical  solutions  of  Eq.  (FI)  shows  that  the  spatial/ temporal  develop¬ 
ment  of  the  phase  distribution  proceeds  gradually  after  an  initial  short  period  of  rapid  changes.  Thus, 

lex )  -  \0„{r  +  1)  -  0„(t)|  «  1,  t  »  1,  (F12) 


ISl  -  y|A0„+l(T)  -  A0„(r)l  «  1.  r  »  1. 


From  Eqs.  (F12)  and  (F5), 


tan  a  =  a  —  0n( r  +  1)  —  9n(r). 
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From  Eqs.  (FI3)  and  (F7), 


sin  S  =s  5  —  —  lAfln+i(r)  —  A 0„(t)], 


From  Eqs.  (FI 3)  and  (F6), 


cos  5=1. 


A~  j  [A 0n+l(r)  +  A 0„(t)]  =  A®b(t). 


Substituting  Eqs.  (F14)  through  (F17)  into  Eq.  (Fll)  we  obtain 


l 0„(t  +  1)  -  0„(t)  1  = 


F  cos  A 0. 


[A9„+| (t)  -  A0„(t)1. 


1 +2FcosAeJl“'"’+1'" 


Using  our  definition  for  A0„(t),  Eq.  (F3),  Eq.  (F18)  may  be  written  as 

[0„(r  +  l)  —  0„(t)1  *=  yC,(r)[^+|(r)  -  20„(t)  +  0b_,(t)),  (F19) 

where 

-i 

c„(r)  -  l+|r'  sec  A«,(t)  .  (F20) 

Equation  (F19)  has  the  form  of  a  1-D  heat  flow  equation  in  which  the  lime  and  space  derivatives  have 
been  replaced  by  finite  differences.  However,  the  solution  of  Eq.  (F19)  is  complicated  by  the  fact  that 
the  "diffusivity"  C„(r)  is  a  function  both  of  space  and  of  lime.  Further  simplification  is  achieved  by 
remembering  that  when  r  »  1  the  quantities  A 0„(r)  approach  their  equilibrium  values, 

lim  A0„(r)  «*  A<£,  (F21) 

T— oo 

where  we  simplify  notation  by  setting  j  =  0  in  Eq.  (19).  From  Eqs.  (F19),  (F20),  and  (F21) 


9„(r  +  1)  -  jC[0„+,(t)  +  +  (1  -  C)0„(t), 


where 


c-  |l  +  -i-  r'  sec  A*J  .  (F23) 

In  the  limit  F  — *  °°  we  see  that  C  —  1,  and  Eq.  (F22)  reduces  as  expected  to  Eq.  (A  10).  More  gen¬ 
erally,  we  solve  Eq.  (F22)  by  the  same  separation-of-variables  method  we  previously  used  to  solve  Eq. 
(A10),  with  the  result 

*Jr)  -  *£  AkZ^,il.  (F24) 

*-i 


where 


m*-CA*  +  (!-C) 


and  where  1 1>„,  Xk  and  Z,*’  are  given  by  Eqs.  (C5),  (C16),  and  (C17).  Again,  from  Eqs.  (F23)  and 
(F25) 

lim  nk  -  Xk,  (F26) 

F— 00 

and  Eq.  (F24)  reduces  to  our  previous  result,  Eq.  (Cl 9). 


41 


RICHARD  A.  STEINBERG 


Assuming  that  the  initial  phase  distribution  is  given  by  Eq.  (C2)  (cf.  also  Fig.  5)  the  expansion 
coefficients  Ak  in  Eq.  (F24)  are  once  again  given  by  Eq.  (C23).  Analogous  to  Eqs.  (C24)  and  (C25)  we 
now  obtain 


«M*)  -  X 

where 


(F27) 


^‘*’(7)  =  lA0n(«>)  -  A<601(-1)*  cot 


1  —  2  C  sin2 


irk 

2M 


(F28) 


Equation  (F27)  is  the  exact  solution  of  the  phase-diffusion  equation,  Eq.  (F22).  However,  Eq.  (F22)  is 
itself  only  an  approximation  to  Eq.  (FI),  valid  in  the  limit  t  »  1.  Thus,  even  if  all  ( M  —  1)  terms 
are  kept  in  our  expansion  for  Eq.  (F27),  the  result  is  still  just  an  approximation.  By  contrast, 

Eq.  (C24)  provides  the  exact  solution  to  Eq.  (Cl). 


Another  distinction  between  Eqs.  (F27)  and  (C24)  is  in  the  eigenvalue  spectrum  that  determines 
the  time-dependence  of  the  space/time  modes,  ili{nkHr).  Previously,  we  found  that  the  A;-mode  was 
degenerate  with  the  ( M  -  Ar)-mode,  i.e.,  that 

(F29) 

Consequently,  it  was  necessary  to  retain  the  highly  oscillatory  ( M  —  l)-mode  in  developing  an  approxi¬ 
mation  to  Eq.  (C24).  However,  we  now  find  that  self-coupling  has  broken  the  degeneracy, 

Im*I  >  <F30) 


and  that  we  need  retain  only  the  slowly  varying  k  “  1  term  in  developing  an  approximation  to  Eq. 
(F27).  Assuming  that  (ir/M)  «  1, 

9„(t)  -  0„(oo)  -  (2A//ir)lA0,,(«>)  -  A^0]  sin  ( nir/M)e  r/T‘,  (F31) 


analogous  to  Eq.  (33).  However,  we  now  have  a  different  value  for  the  settling  time, 

1  +  y  r'  sec  A<J. 


(F32) 


where  rs  —  A/2/ 5  is  the  settling  time  for  modules  without  self-coupling,  Eq.  (34).  In  general,  when  F 
>  0,  stable  equilibrium  phase  distributions  must  have  0  <  A</>  <  rr/2;  when  F  <  0,  we  have  w/2  < 
A <t>  <  it.  It  then  follows  from  Eq.  (F32)  that 


Ts  ^  t5, 

always;  i.e.,  modules  without  self-coupling  settle  faster  than  modules  with  self-coupling. 


(F33) 


Instead  of  Eq.  (35)  we  now  obtain 

A®, (r)  -  A0b(<»)  -  2  (A0„(<»)  -  A$01  cos  ( nn/M )  e  r/T‘,  (F34) 

where  Ts  is  once  again  given  by  Eq.  (F32).  Comparing  the  n-dependence  of  Eqs.  (35)  and  (F34)  we 
see  that  Eq.  (F34)  is  missing  the  highly  oscillatory  factor  p  that  represents  the  contribution  of 
to  Eq.  (35). 


We  have  compared  the  one-mode  approximation,  Eq.  (F34),  with  some  exact  computer  solutions 
to  Eq.  (FI).  In  performing  these  calculations  we  assume  that 

(A/.  A0O.  aMm  -  (9,  60°,  45°)  (F35) 

as  for  Eq.  (36)  and  Fig.  8.  Figure  FI  presents  a  comparison  between  the  analytic  and  numerically 
derived  phase  distributions  as  a  function  of  the  module  index  n  for  two  particular  values  of  time 
(7  —  70  and  7  -  100). 
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Compared  with  Fig.  8  we  note  from  Fig.  FI  that  the  phase  distribution  for  self-coupled  modules 
evolves  much  more  slowly  with  time  and  is  spatially  smoother  than  the  phase  distribution  of  modules 
that  lack  self-coupling. 
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